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Executive  Summery  i  The  characteristics  of  the  ambient  noise  field  in 
the  sea  are  important  to  the  performance  prediction,  detection  capabilities 
and  optimum  use  of  sonar  system;.  Rsrngniring  the  aim  of  ambient 
noise  stndies  at  SAOlrANTCBlf^s  to  improve  knowledge  and  prediction  of 
the  ambient  noise  Held  in  areas  of  importance.  These  studies  involve  the 
development  and  validation  of  models  for  the  generation  and  distribution  of 
the  noise. 


This  report  presents  a  model  for  the  generation  of  high-frequency,  wind- 
related  noise.  The  model  is  based  on  the  assumption  that  the  noise  is 
generated  by  air  bubbles  entrained  in  whitecaps  (breaking  waves).  The 
characteristics  of  the  noise  predicted  by  the  model  compare  favourably  with 
published  observations,  regarding  directionality,  level  and  spectral  form. 

The  development  of  the  model  involves  a  statistical  approach  in  terms  of 
the  individual  bubble*  rather  than  a  continuum  approach  in  terms  of  an 
alr/water  mixture. 

This  work  concludes  the  theoretical  modelling  of  high-frequency,  wind- 
related  noise ,csndmeted  at  the  CentW^The  results  contribute  to  progress  in 
identifying  the  sources  of  wind-related  noise  in  the  sea.  A  later  report  wiii 
present  the  results  of  experimental  observations  of  the  noise  from  whitecaps 

at  sea. 
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A  model  for  the  generation  of  noise  by 
bubbles  formed  in  breaking  waves 

R.M,  Heitmeyer  and  R,D.  Hollett 


Abstract:  A  general  model  is  presented  that  describes  the  high-frequency 
noise  cross-spectrum  generated  by  the  vibrations  of  bubbles  formed  in  break¬ 
ing  waves.  The  bubble  vibrations  are  described  s  linear  with  the  excitation 
occurring  at  the  moment  the  bubble  is  formed.  Two  simplified  models  ate 
obtained  from  the  general  model  under  idealized  propagation  assumptions. 
The  simplified  models  can  be  viewed  as  extensions  of  the  Kupennan  and 
Ingenito  and  the  Cron  and  Sherman  models  in  that  the  source  mechanism  is 
included.  To  obtain  specific  results  we  postulate  that  the  bubble  vibrations 
are  excited  either  by  a  rapid  pressure  change  at  formation  or  by  an  initial 
rate  of  change  of  the  bubble  volume.  The  noise  intensity  spectra  for  these 
two  excitation  mechanisms  are  evaluated  over  the  frequency  range  0.66-10 
kHz  using  one  of  the  simplified  models  with  the  parameters  estimated  for  a 
particular  wind  speed.  Although  different,  both  spectra  have  levels  compa¬ 
rable  to  the  Wene  spectrum  for  the  same  wind  speed. 


Keywords:  breaking  waves  o  bubble  vibrations  o  wind-related  noise 


SaCLANTCEN  SR- 145 


Contents 

1.  Introduction .  1 

2.  The  general  model .  2 

3.  The  simplified  models .  8 

4.  The  excitation  postulates . 13 

5.  Parameter  estimation  and  model  results . 14 

6.  Summary  and  conclusion  . 20 

References  . 21 

Appendix  A  -  The  breaking-wave  and  the  bubble  occurrences  .  .  23 

Appendix  B  -  Derivation  of  the  noise  cross-spectrum . 26 

Appendix  C  -  Derivation  of  the  breaking-wave  cross-spectrum  .  ,  30 

Appendix  D  -  The  high-frequency  approximation . 34 

Appendix  E  -  The  model  equations  for  the  examples  of  Sect.  5  ..  38 


SACLANTCEN  SR-145 


1.  Introduction 

In  his  classical  paper,  Wenz  [1]  speculated  that  two  likely  sources  of  high-frequency, 
wind-related  noise  are  spray  and  bubbles.  Spray  is  generated  by  the  action  of  the 
wind  on  the  crests  of  waves  and  by  wave  breaking.  Bubbles  are  generated  in  large 
numbers  by  the  entrainment  of  air  in  breaking  waves.  Air  entrainment  following 
the  impact  of  spray  also  generates  bubbles.  A  number  of  authors  have  developed 
models  for  the  high-frequency  noise  generated  by  spray  impact  or  bubble  vibra¬ 
tions.  The  spray  impact  model  of  Wilson  [2]  is  an  empirical  model  based  on  Franz’s 
measurements  [3]  of  the  sound  generated  by  the  impact  of  sprays  in  a  water  tank. 
To  rigorously  quantify  the  relative  contribution  of  spray  impact  noise,  further  work 
must  be  done  on  the  amount  of  spray  generated  and  the  characteristics  of  that  spray 
(droplet  size,  impact  speed  and  angle).  Kerman  [4]  has  developed  a  model  for  the 
noise  generated  by  non-linear  vibrations  of  the  bubbles  in  breaking  waves.  The  as¬ 
sumption  is  made  that  the  turbulence  in  breaking  waves  is  sufficiently  strong  to  force 
non-linear  vibrations  of  the  entrained  bubbles,  Other  models  that  are  also  based  on 
non-linear  motion  of  the  bubble  wall  are  those  of  Furduev  [5]  and  Shang  and  Ander 
son  [6].  Whether  or  not  non-linear  motion  of  the  bubble  wall  really  occurs  remains 
to  be  determined.  Crowther  [7]  has  developed  a  semi-empirical  model  in  which  the 
noise  is  generated  by  linear  vibrations  of  the  bubbles  excited  at  formation.  In  his 
model  the  source  level  is  determined  empirically  from  the  noise  spectra  of  bubbles 
formed  followivg  spray  impact  in  a  tank. 

In  this  report  we  present  a  model  for  the  noise  generated  by  the  linear  vibrations  of 
the  bubbles  formed  in  breaking  waves  where  the  vibrations  are  assumed  to  be  excited 
at  the  moment  the  bubble  is  formed.  General  expressions  are  presented  for  the  cross¬ 
spectra  of  the  noise  produced  by  an  individual  breaking  wave  and  the  noise  produced 
by  the  aggregate  of  the  breaking  waves.  The  spectra  for  a  single  hydrophone  are 
obtained  as  special  cases  of  the  cross-spectra.  We  then  present  two  simplified  models 
for  the  cross-spectra  obtained  under  idealized  propagation  assumptions.  Both  the 
general  and  the  simplified  models  are  applicable  to  an  arbitrary  bubble  excitation 
mechanism.  To  obtain  specific  results  we  postulate  that  the  bubble  vibrations  are 
excited  either  by  a  rapid  change  in  the  pressure  on  the  bubble  wall  at  the  moment  of 
formation  or  by  an  initial  rate  of  change  of  the  bubble  volume.  The  noise  intensity 
spectra  for  these  two  excitation  mechanisms  are  evaluated  over  the  frequency  range 
0.66-10  kHz  using  one  of  the  simplified  models  and  ihe  results  are  compared  to  a 
Wenz  spectrum. 
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2.  The  general  model 


The  sound  produced  by  the  linear  vibrations  of  a  bubble  is  well  understood.  The 
pressure  waveform  is  determined  from  the  second  time  derivative  of  the  volume 
vibrations  which  in  turn  are  determined  as  the  solution  of  a  second-order  linear 
differential  equation  with  constant  coefficients  [8,9].  Accordingly,  the  spectrum  of 
the  pressure  observed  at  some  reference  distance  *o  can  be  written  as  the  product 

a(*)=ffMF(u»),  (1) 


where  H(u>)  is  the  pressure  transfer  spectrum  and  F(u>)  is  the  pressure  excitation 
spectrum.  The  pressure  transfer  spectrum  is  equal  to  the  bubble  pressure  spectrum 
when  the  bubble  is  excited  by  an  impulse  excitation  of  unit  strength.  This  spectrum 


is  given  by 


-{It/roW1 

(u>J  -  )  4-  tiu i0w  ' 


(2) 


where  R  is  the  equilibrium  radius  of  the  bubble,  u)0  ~  2 t/0,  is  the  natural  angular 
frequency  of  the  monopole  bubble  vibrations,  and  5  is  the  damping  constant.  For 
bubble  radii  of  interest  here,  u>0  and  S  are  well  approximated  [9,10]  by 


oj0  =  2ir(3.3  m  s  1 )/ R,  (3a) 

6  =  0.014+  1.8  x  10  4u4/2.  (3b) 


The  pressure  transfer  spectrum  for  the  radii  of  interest  is  dominated  by  a  large  peak 
that  occurs  for  u>  ss  w0,  with  a  value  approximately  given  by  R/Sx0.  The  pressure 
excitation  spectrum  describes  both  the  applied  pressure  on  the  bubble  wall  and  the 
departure  from  the  equilibrium  state  at  the  time  at  which  the  bubble  is  caused  to 
vibrate.  This  spectrum  is  given  by 

F{w)  =  ( p/4rrR )  (t'wAV(O)  +  V(0))  -  AP( u>),  (4) 

where  AV(0)  =  V(0)  -  Vo,  is  the  initial  volume  displacement,  Vo  =  is  the 

equilibrium  volume  of  the  bubble,  V(0)  is  the  initial  volume  rate,  and  AP(w)  is 
the  spectrum  of  the  pressure  excess  acting  on  the  bubble  wall.  The  initial  volume 
displacement  is  equivalent  to  an  excitation  by  a  doublet  (the  time  derivative  of  an 
impulse),  and  the  initial  volume  rate  is  equivalent  to  an  impulse  excitation.  Thus,  the 
first  term  has  a  significant  spectral  component  at  the  natural  bubble  frequency  and 
hence  gives  rise  to  high-frequency  sound.  This  term  is  non-zero  whenever  a  bubble  is 
formed  in  a  non-equilibrium  state,  The  pressure-excess  term  has  significant  energy 
at  the  natural  bubble  frequency  if,  at  the  moment  of  formation,  the  excess  pressure 
waveform  changes  significantly  over  a  time  period  of  the  order  of  2x/wo  Whether 
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the  bubble  is  excited  by  formation  in  a  non-equilibrium  state  or  by  a  rapid  change 
in  excess  pressure  at  formation,  the  resulting  excess  energy  is,  in  part,  emitted  in 
the  form  of  a  damped  pressure  wave.  The  spectrum  of  the  pressure  wave  observed 
at  any  point  is  determined  as  the  product  of  the  spectrum  4irj;o^>b(w)  and  the 
appropriate  Green  function.  In  what  follows  we  omit  the  reference  distance  *0  with 
the  understanding  that  unit  distance  is  implied. 

A  physical  picture  for  the  formation  of  the  bubbles  that  produce  the  sound  is  pro¬ 
vided  by  the  description  of  a  spilling  breaker  by  Longuet- Higgins  and  Turner  [11], 
According  to  these  authors,  the  whitecap  formed  on  a  spilling  breaker  can  be  re¬ 
garded  as  a  turbulent  air/water  mixture  that  is  both  accelerated  down  the  forward 
slope  of  the  wave  by  gravity  and  retarded  by  the  entrainment  of  upslope  momentum 
from  the  flow  in  the  wave  below.  This  air/ water  mixture  is  lighter  than  the  water  in 
the  wave  below  owing  to  the  entrainment  of  air  bubbles  and  remains  distinct  from 
the  rest  of  the  wave.  Bubbles  are  generated  at  the  front  of  the  whitecap  by  the 
over-running  of  a  layer  of  air  as  the  front  advances  and  they  are  generated  along 
the  surface  of  the  whitecap  by  the  trapping  of  air  under  the  turbulent  eddies  which 
break  out  of  the  surface  of  the  whitecap.  We  consider  only  those  bubbles  that  are 
generated  at  the  whitecap  front  -  not  those  that  are  generated  along  the  surface  of 
the  whitecap  itself. 

The  model  for  the  cross  spectrum  of  the  noise  generated  by  a  breaking  wave  is 
based  on  two  fundamental  assumptions.  Firstly,  neither  the  sound  produced  by  an 
individual  bubble  nor  the  propagation  of  that,  sound  to  the  hydrophones  is  influenced 
by  the  presence  of  the  other  bubbles  in  the  breaking  wave.  Stated  in  other  terms, 
the  bubbles  in  the  breaking  wave  do  not  vibrate  collectively  to  produce  the  sound 
and  the  propagation  of  the  sound  from  an  individual  bubble  is  not  absorbed  or 
scattered  by  the  other  hubbies  in  the  breaking  wave.  Secondly,  the  number  of  bubbles 
generated  in  the  breaking  wave  and  the  positions  and  times  at  which  those  bubbles 
are  generated  are  described  by  a  Poisson  process.  Furthermore  the  equilibrium 
radii  of  those  bubbles  are  statistically  independent  with  a  probability  density  that 
depends  oidy  on  the  bubble  occurrence  position. 

The  setting  for  the  model  is  illustrated  in  Fig.  1.  The  breaking  wave  occurs  at  posi¬ 
tion  r  and  time  r  and  the  sound  produced  by  that  wave  is  observed  on  hydrophones 
located  at  positions  and  zj  The  bubbles  are  generated  in  a  region  located  at 
the  front  of  the  whitecap  that  advances  with  the  whitecap  front  down  the  leading 
face  of  the  breaking  wave.  The  total  volume  swept  out  by  the  bubble  generation 
region  during  the  lifetime  of  the  breaking  wave  is  the  source  region  V .  The  total 
number  of  bubbles  generated  in  the  source  region  is  N\,  and  the  positions  and  times 
at  which  those  bubbles  are  generated  are  {(y*,  v^),  k  —  1, . . . ,  N t,},  where  is  mea¬ 
sured  relative  to  the  breaking-wave  occurrence  time  r.  The  damped  pressure  wave 
emitted  by  each  bubble  is  of  sufficiently  short  duration  that  the  effect  of  buoyancy 
and  transport  of  the  bubble  in  the  turbulence  is  neglected.  By  the  first  assumption 
the  spectrum  of  the  pressure  observed  on  each  hydrophone  from  all  of  the  bubbles 
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Fig.  1.  The  physical  setting  fot  the  general  model. 


generated  in  the  breaking  wave  is  the  superposition  of  the  received  pressure  spectra 
from  each  of  the  bubbles.  Thus  using  the  linear  theory  to  describe  the  received 
pressure  spectrum  from  each  bubble,  the  spectrum  of  the  total  pressure  observed  on 
each  hydrophone  can  be  written  in  the  form 


ATfc 

Pw(u>;r,Zj )  =  V  Ph(u-yk,Rk)G[u,r  +  yfc, 

fc=i 


(5) 


where  P\,(ur,y,  ?)  is  the  pressure  spectrum  for  a  bubble  generated  at  position  y 
with  equilibrium  radius  R  and  G(u*,r  -f  y,ar,)/4jr  is  the  Green  function  describing 
the  propagation  from  the  bubble  position  r  +  y  to  the  hydrophone  positions  zi  and 
Zj .  For  fixed  values  of  y,  v  and  R,  the  Green  function  and  the  pressure  transfer  spec¬ 
trum  H(w,y,R)  are  assumed  to  be  deterministic;  the  pressure  excitation  spectrum 
F(w;y,  R)  is,  in  general,  stochastic  with  cross-spectrum 

Gr(u,;y,y\ £,£')=  Eh{F‘(u,;y,  R)F(w,y' ,  R‘)\  (6a) 

and  energy  spectrum 

Sr(u\y,R)=  £b[|F(w;y,il)|*],  (6b) 

where  £/’b[  ]  is  the  expected  value  operator  for  fixed  values  of  the  bubble  occurrences 
and  the  bubble  radii. 

The  Poisson  process  describes  the  generation  of  bubbles  by  the  ensemble  of  breaking 
waves  that  occur  at  the  position  r.  This  process  is  specified  by  a  bubble  generation 
rate  n(y,v)  that  represents  the  mean  number  of  bubbles  that  are  generated  per 
unit  volume  per  unit  time  in  the  bubble  generation  region  during  the  breaking- wave 
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lifetime.  The  mean  number  generated  per  unit  vohune  at  the  point  y  in  the  source 
region  V  is  given  by  the  bubble  generation  rate 

-  /  My.  ")<*<'.  (7») 

J  0 

and  the  mean  number  of  bubbles  is 

Mb  =  /  My)d3y,  (7b) 

Jv 

where  d3y  is  the  volume  element.  Finally  the  equilibrium  radii  of  the  bubbles 
generated  at  the  position  y  are  described  by  the  bubble-radius  probability  density 
p(R  |  y).  For  reference  purposes,  the  properties  of  Poisson  processes  relevant  to  this 
report  are  summarized  in  Appendix  A. 


The  breaking-wave  cross-spectrum  is  the  mean  cioss-spectrum  of  the  pressure  ob¬ 
served  at  Zi  and  z 2  for  all  breaking  waves  that  occur  at  the  position  r.  In  Appendix  C 
it  is  shown  that  this  cross-spectrum  can  be  written  as  the  stun  of  two  cross-spectra. 
The  first  is  the  cross-spectrum  that  results  if  the  received  pressure  contributions 
from  the  individual  bubbles  were  statistically  uncorrelated  with  zero  mean.  The 
second  cross-spectrum  includes  the  effects  of  both  the  non-zero  mean  and  the  cor¬ 
relation  in  the  received  pressure  contributions  that  results  from  the  correlation  in 
the  bubble  excitation  spectra  (see  (6a))  In  Appendix  D  we  derive  an  upper  bound 
on  the  energy  contribution  of  the  second  spectrum  relative  to  that  of  the  first.  Dy 
evaluating  this  upper  bound  we  show  that  the  first  spectrum  makes  the  dominant 
contribution  at  the  high  frequencies  of  interest  in  this  report.  On  the  basis  of  this 
result  we  henceforth  neglect  the  contribution  of  the  second  cross-spectrum.  With 
this  approximation,  the  breaking- wave  cross-spectrum  can  be  written  in  the  form 


Cw(w;r,zi,z2)  =  MbSm  b  (w)Cp(w;  r,  zt)  z2), 


(8a) 


where  S'm.biw)  is  the  mean  bubble  spectrum  and  Cp(u>;r,zi,z2)  is  the  propagation 
cross-spectrum.  The  mean  bubble  spectrum  is  given  by 


Sm.b.(w)  =  Mb  1  [  p(y)Sth.(w,y)diy, 
Jv 


where 


St.b  (w;y)  =  f  Sb{v,y ,  R)p(R  |  y)dR 
Jo 

is  the  radius- averaged  bubble  spectrum  and 

Sb(^y1R)  =  \H(w,y,R)\iSF(u,-,y,R) 
is  the  single-bubble  spectrum.  The  propagation  cross-spectrum  is  given  by 

Cp(w;r,z,,z2)  =  /  Sr.b.(u;;y)6’*(w,r  +  y,zi)G(w,r  +  y,z2)d3y, 
Jv 


(8b) 

(8c) 

(8d) 

(8e) 
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where 


St.  b.(w;y) 


/<(y)S,.b.(a>;y) 


,b.(w;y)d3y 


(8f) 


is  the  normalized  radius-averaged  spectrum.  The  complete  equation  for  the  breaking- 
wave  cross-spectrum  car  be  found  in  Appendix  C. 


The  cross-spectrum  of  (8a)  can  be  interpreted  as  follows.  The  propagation  cross¬ 
spectrum  Cp(u;;r,  zi,z2)  is  the  cross-spectrum  of  the  pressure  from  a  volume  dis¬ 
tribution  of  uncorrelated,  zero-mean  monopole  sources  where  the  energy  emitted 
by  each  source  is  described  by  the  normalized  radius- averaged  spectrum  5r,t>.(w;y). 
By  virtue  of  (8f)  the  energy  spectrum  level  for  the  volume  distribution  is  normal¬ 
ized  to  unity.  The  single-bubble  spectrum  S\,(w,y,R)  describes  the  energy  emitted 
by  a  bubble  that  is  generated  at  position  y  with  an  equilibrium  radius  R.  The 
radius-averaged  spectrum  St.b.(a>;y)  is  the  average  of  the  single-bubble  spectrum 
with  respect  to  the  radius  probability  density  p(R  |  y);  i.e.  it  describes  the  average 
energy  emitted  by  a  single  bubble  located  at  the  position  y  as  the  radius  of  that 
bubble  ranges  through  all  possible  values.  The  quantity  Mh~1/i(y)  is  the  probability 
density  on  the  bubble  occurrence  positions  (see  Appendix  A).  The  mean  bubble 
spectrum  Sm.b.(w)  is  the  average  of  Sr.b.(^;y)  with  respect  to  this  probability  den¬ 
sity.  Thus  the  mean  bubble  spectrum  describes  the  average  energy  emitted  by  a 
single  bubble  located  anywhere  in  the  source  region  V  as  both  the  radius  and  the 
position  of  that  bubble  range  through  all  possible  values.  Finally  since  M\,  is  the 
mean  number  of  bubbles  generated  in  the  source  region,  the  spectrum  Mb5m. b.(w) 
describes  the  average  energy  emitted  by  all  of  the  bubbles  in  a  breaking  wave.  With 
these  interpretations  (8a)  indicates  that  the  breaking-wave  cross-spectrum  is  deter¬ 
mined  as  the  product  of  a  spectrum  that  describes  the  average  energy  of  the  sound 
emitted  by  a  breaking  wave  and  a  cross-spectrum  that  describes  the  propagation  of 
that  sound  from  the  position  of  the  breaking  wave  to  the  hydrophone  positions. 


The  noise  cross-spectrum  is  defined  as  the  fourier  transform  of  the  cross-correlation 
function  of  the  pressure  due  to  the  aggregate  of  the  breaking  waves  on  the  ocean 
surface.  The  expression  for  the  noise  cross-spectrum  is  obtained  under  two  assump¬ 
tions:  firstly,  that  the  pressure  waveforms  received  at  the  hydrophones  from  dif¬ 
ferent  breaking  waves  are  statistically  independent  and  identically  distributed;  and 
secondly,  that  the  occurrence  positions  and  occurrence  times  of  the  breaking  waves 
{(r*,  T*)}  are  described  by  a  second  Poisson  process  that  is  independent  of  the  bub¬ 
ble  generation  process.  This  process  is  determined  by  a  breaking- wave  occurrence 
rate  A(r)  that  represents  the  mean  number  of  breaking-wave  occurrences  per  unit 
area  of  the  ocean  surface  per  unit  time  at  the  position  r.  From  these  assumptions 
it  is  shown  in  Appendix  B  that  the  noise  cross- spectrum  is  given  by 

C(u/;zi,z2)  =  [  A(r)Cw(w;r,Zi,z2)  dJr,  (9) 

J  A, 
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where  d2r  is  the  surface  area  element  and  A ,  is  the  total  area  of  the  ocean  sur¬ 
face  over  which  the  breaking  waves  occur.  This  equation  indicates  that  the  noise 
cross-spectrum  is  obtained  as  the  integral  over  the  ocean  surface  of  the  mean  cross- 
spectrum  of  the  pressure  from  a  single  breaking  wave  weighted  by  the  breaking-wave 
occurrence  rate.  Note  that  this  result  is  the  same  as  would  be  obtained  if  the  field 
of  individual  breaking  waves  were  viewed  as  a  continuum  on  the  ocean  surface.  It  is 
emphasized,  however,  that  (9)  is  a  consequence  of  the  Poisson  assumption  and  does 
not  result  from  a  continuum  approximation.  Such  an  approximation  would  not  be 
plausible  in  the  light  of  the  sparsity  with  which  breaking  waves  occur  on  the  ocean 
surface. 

For  the  high  frequencies  of  interest  in  this  report  the  noise  cross-spectrum  is  well 
approximated  by  using  (8)  to  substitute  for  the  breaking-wave  cross-spectrum  in 
(9).  The  resulting  noise  cross-spectrum  can  be  written  in  the  form 


C(ur,zl,z2)  =  MhSm.b.(w)Ci.p.(<ir,Zi,z3),  (10a) 

where 

<7i.p.(w;zi,z2)  =  (  A(r)Cp(w;r,z,,z2)d2r  (10b) 

J  A9 

is  referred  to  here  as  the  integrated  propagation  cross-spectrum.  To  interpret  these 
equations  we  note  that  the  noise  cross-spectrum  can  be  expressed  as  the  integral  over 
the  surface  region  of  the  product  of  the  spectrum  A(r)Jl/|,Sm.b  (u>)  and  the  propaga¬ 
tion  cross-spectrum  C’p(u.';r,Zi,z2).  As  noted  previously  the  spectrum  MbSm.b.(n>) 
is  the  mean  energy  spectrum  for  the  breaking  waves  that  occur  at  the  position  r, 
and  A(r)  is  the  rate  per  unit  surface  area  per  unit  time  at  which  these  breaking 
waves  occur.  Jt  follows  that  the  product  A(r)Mb5m.b.(w)  is  the  mean  intensity 
source  spectrum  per  unit  area  for  the  succession  of  breaking  waves  that  occur  at  r. 
Consequently,  the  product  of  this  source  spectrum  with  the  propagation  spectrum 
Cp(w;r,  zi,z2)  is  the  noise  cross-spectrum  per  unit  area  for  the  breaking  waves  that 
occur  at  r.  Thus  (10a)  and  (10b)  state  that  the  noise  cross-spectrum  for  all  breaking 
waves  is  obtained  as  the  integral  of  the  noise  cross-spectrum  for  the  breaking  waves 
that  occur  at  each  point  on  the  surface  region. 
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3.  The  simplified  models 


In  the  noise  model  the  propagation  of  the  sound  to  the  hydrophone  positions  is 
described  by  the  two  cross-spectra  Cp(w;  r,  zj ,  z2 )  and  Ci.p.(u>;  zi ,  z2).  The  propaga¬ 
tion  cross-spectrum  Cp(u/;r, zt,z2)  describes  the  sound  propagation  from  all  points 
in  the  source  region  of  a  breaking  wave  located  at  the  position  r;  the  integrated 
cross-spectrum  C j.p.(u>;  zi .  z2)  describes  the  sound  propagation  from  all  breaking- 
wave  positions.  A  considerable  simplification  in  the  noise  model  results  if  the  acous¬ 
tic  propagation  is  appoximated  as  that  determined  for  a  source  region  lying  below  a 
flat,  pressure-release  surface.  With  this  approximation  the  integrated  cross- spectrum 
can  be  obtained  from  ah  extension  of  the  noise-propagation  model  of  Kuperman  and 
Ingenito  [12],  provided  that  the  source  region  has  a  horizontal  cross-section  that  is 
independent  of  the  vertical  coordinate,  and  that  the  medium  is  horizontally  stratified 
(i.e.  the  sound  speed  and  bottom  characteristics  are  independent  of  the  horizontal 
coordinates).  If,  in  addition,  the  medium  can  be  viewed  as  homogeneous  and  semi¬ 
infinite  (constant  sound  speed  and  no  bottom)  the  integrated  cross-spectrum  can  be 
obtained  from  an  extension  of  the  Cron  and  Sherman  model  [13].  In  this  section  we 
describe  the  simplified  models  that  result  from  these  approximations. 

The  setting  for  the  simplified  models  is  illustrated  in  Fig.  2.  We  take  the  surface 
region  As  to  be  of  infinite  extent  and  assume  that  the  breaking-wave  occurrence 
rate  A(r)  is  independent  of  r.  The  source  region  extends  down  from  the  surface  to  a 
depth  D  with  a  horizontal  cross-section  .4'  that  is  independent  of  depth.  Within  the 
source  region,  the  bubble  position  vector  y  is  represented  as  the  sum  of  a  horizontal 
vector  y'  and  a  vertical  vector  y"  with  magnitude  y"  £  [0.  D).  The  two  hydrophones 
are  located  at  positions  zi  and  z2.  In  the  ideal  medium  case  (homogeneous,  semi¬ 
infinite)  the  hydrophones  are  assumed  to  be  located  far  enough  below  the  ocean 
surface  that  the  incident  acoustic,  field  can  be  viewed  as  planewave.  For  this  case, 
the  vertical  vector  z,  of  magnitude  z,  points  to  the  midpoint  of  the  line  connecting 
the  hydrophones.  The  positions  of  the  hydrophones  relative  to  the  midpoint  are 
determined  by  the  hydrophone  diiference  vector  z'  -  z2  --  z\ ,  that  has  magnitude 
z'  and  lies  at  an  angle  7  relative  to  the  horizontal  plane.  The  planewave  field  is 
described  by  the  wavenumber  vector  k  that  has  magnitude  k  —  u>/c,  and  points  in 
the  direction  of  z  -  r  at  an  angle  9  relative  to  the  vertical. 

In  the  stratified  medium  case,  the  noise  cross-spectrum  is  obtained  from  the  gen¬ 
eral  model  using  the  integrated  cross-spectrum  that  results  from  the  simplifying 
assumptions.  To  obtain  the  integrated  cross-spectrum,  we  use  the  propagation 
cross-spectrum  of  (8e)  in  (10b),  interchange  the  order  of  the  surface  and  the  vol¬ 
ume  integrations,  and  express  the  volume  integral  a3  an  integral  over  the  horizontal 
cross-sectional  area  followed  by  an  integral  over  the  depth  interval.  (This  last  step 
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Fig.  2.  Setting  for  the  simplified  models. 


is  permitted  by  the  assumption  on  the  form  of  the  source  region.)  The  result  is 

Ci.p.(u>; *!,**)  =A  f  I  5,.b.(u>;y'  -i-  y") 

Jo  J  A' 

x  f  G'(w,r  +  y'  +  y",*i)G(u),r  +  y'  +  y", *2)d3rdsy.  (11) 
Ja, 

The  final  expression  for  the  integrated  cross-spectrum  is  obtained  by  making  the 
change  of  variable  r'  =  r  +  y'  in  the  surface  integral  of  (11)  and  using  the  assumption 
that  At  tas  infinite  extent  to  eliminate  the  dependence  on  the  horizontal  vector  y'. 
The  noise  cross- spectrum  is  then  obtained  by  substituting  in  (10a).  The  result  can 
be  written  in  the  form 


C{  w;zi,z2)  =  AMbSm.b.MCftwjzi.z:), 


(12a) 


where 


[D  ■ 

Cl(w,zuzi)=  I  S,.k.(w\y")C'{w\y",zuz2)dy"  (12b) 

Jo 


and 


^r.b.(w;y")=  /  •S,rb.(u>;y'  +  y")d2y', 

J  A' 

C"(w;y".*i.*:)  =  /  G*(^,r  +  y",z1)G(w,r  + y",*2)d:i 
J  A, 


(12c) 

(12d) 
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Note  that  C|(tu;zi,*2)  is  the  integrated  cross-spectrum  that  results  for  a  breaking- 
wave  occurrence  rate  of  unity. 

The  simplified  model  of  Eqs.  ( 12)  determines  the  noise  cross- spectrum  as  the  product 
of  the  source  spectrum  AMbSm.b.(w)  and  the  integrated  cross-spectrum  C{(  u>;  z\ ,  z2). 
According  to  (12b),  the  integrated  cross-spectrum  is  obtained  as  a  depth  average 
of  the  cross-spectrum  C'(u>;  y",  Z\ ,  z2)  weighted  by  the  depth- dependent  spectrum 
b.(w;  y").  The  cross-spectrum  C'(u>;  y",  zi ,  z2)  is  equivalent  to  that  obtained  in 
the  Kuperrnan  and  Ingenito  model  for  a  sheet  of  spatially-uncorrelated  monopole 
sources  located  at  a  depth  y"  below  the  surface.  (To  see  this,  make  the  change  of 
variables  p  —  r'  -  r"  in  Eq.  (7)  of  [12],  and  take  the  sources  to  be  uncorrelated  by 
using  Eq.  (17)  of  [12],  to  obtain  (12d).)  Thus  the  above  simplified  model  provides 
both  an  expression  for  the  source  spectrum  of  the  Kuperrnan  and  Ingenito  model 
and  a  rule  for  combining  the  ‘source- sheet’  cross-spectra  determined  by  their  model 
to  produce  the  total  noise  cross-spectrum.  Examples  showing  the  response  of  an 
array  determined  from  the  ‘source-sheet’  cross-spectrum  of  Kuperrnan  and  Ingenito 
for  different  acoustic  environments  can  be  found  in  Hamson  [14]. 

In  the  ideal  medium  case  the  resulting  simplified  model  provides  expressions  for  both 
the  breaking-wave  cross-spectrum  and  the  noise  cross- spectrum.  These  expressions 
are  obtained  from  the  general  model  with  the  propagation  cross-spectrum  and  the 
integrated  cross-spectrum  determined  from  the  Green  function  for  the  ideal  medium. 
Using  the  method  of  images  the  Green  function  can  be  obtained  by  replacing  the 
point  source  in  the  semi-infinite  medium  by  a  dipole  source  in  an  infinite  medium 
where  the  dipole  axis  is  perpendicular  to  the  surface  and  the  length  of  the  dipole  is 
twice  the  bubble  depth.  It  follows  immediately  that  the  function  G(w,r+  y,z>)  is 
given  by 

-,~a|i-r|/2 

G(w.r  +  y  ,Zj)  =  2»  sm(k"y")e~,k‘>  — ; - — ,  (13) 

I*  ~  r| 

where  a  is  a  frequency-dependent  absorption  coefficient  and  k’  is  the  vertical  com¬ 
ponent  of  k.  The  breaking- wave  cross-spectrum  is  obtained  by  first  using  (13)  in 
(8e)  to  determine  the  propagation  cross-spectrum  and  then  substituting  in  (8a).  The 
result  can  be  expressed  in  the  form 


Cw(u>;r,zi,z2)  =  MbSmb.{w)R(w,6)T{u>,  |z  - 

r|)e-<k 

(14a) 

where 

R(u>,0)  =  f  5r,b.(w;  y")  (2  sin  (fr  cos(0)y" 
Jo 

))’  dy" 

(14b) 

is  a  frequency-dependent  radiation  pattern  and 

-a|»-r| 

(14c) 
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is  the  transmission  loss  from  rtos.  This  cross-spectrum  is  the  same  as  that  obtained 
for  a  point  source  that  has  an  energy  source  spectrum  AfbSmb.(u>)  and  a  radiation 
pattern  R(u,Q)  and  that  radiates  sound  into  an  infinite  homogeneous  medium. 

The  noise  cross-spectrum  is  obtained  with  the  integrated  cross-spectrum  determined 
from  the  propagation  cross-spectrum  for  the  ideal  medium.  The  integral  in  (10b) 
is  first  expressed  in  polar  coordinates  to  evaluate  the  azimuth  integral  and  then 
converted  to  an  integral  with  respect  to  $  by  making  the  change  of  variables  |r|/z  = 
tan(0).  The  result  is  then  combined  with  (10a)  to  obtain 

<?(<*>;  2Si,*i)  =  AMb.9n,.b.(u>)<r’.(u>;2,z',7),  tl5a) 

where  the  integrated  cross-spectrum  C'j'fw;  2,  2', 7)  is  given  by 

»*/2 

C|'(u>; 2, r',7)  =  2ir  / 

Jo 

X  J0  (£2'  cos(7)sin(0))tan(0)d0.  (15b) 

This  noise  cross-spectrum  is  the  same  as  that  obtained  by  replacing  the  ocean  sur¬ 
face  with  a  sheet  of  uncorrelated  point  sources  that  radiates  sound  into  an  infinite 
homogeneous  medium  where  each  source  has  an  intensity  source  spectrum  per  unit 
area  AAfb5mb.(u>)  and  a  radiation  pattern  R(u>,0). 

The  simplified  model  of  Eqs.  (15)  is  related  to  the  Cron  and  Sherman  model  through 
the  real  part  of  the  normalized  integrated  cross-spectrum  plw;  2, 2',  ■) )  -  Re  (C((w; 
2,2',  7)/Cj'(w;  2,0,7).  Using  Eq.  (15b)  to  calculate  p(ui;  z,  2', 7)  and  comparing 

the  result  with  Eq.  (14)  of  (13),  it  is  seen  that  for  a  lossless  medium  (a  =  0), 

the  two  cross-spectra  are  identical  if  the  radiation  pattern  R(u>,0)  is  replaced  by 
cos2,n(0).  Thus  this  simplified  model  provides  not  only  a  source  spectrum  for  the 
noise-propagation  model  of  Cron  and  Sherman  but  also  an  expression  for  the  radi¬ 
ation  pattern  in  their  model. 

Further  insight  can  be  obtained  by  comparing  the  angular  dependence  of  the  radi¬ 
ation  pattern  R(t*>,0)  at  extreme  frequencies  with  the  assumed  angular  dependence 
in  the  radiation  pattern  of  Cron  and  Sherman.  To  this  end  it  is  noted  that  the  radi¬ 
ation  pattern  iJ(u>,0)  is  obtained  as  the  depth  integral  of  the  radiation  pattern  for 
a  vertical  dipole  weighted  by  the  depth-dependent  spectrum  5,,b.(w;y").  For  low 
frequencies  the  dipole  radiation  pattern  is  approximately  that  of  an  infinitesimal 
dipole  for  all  depths  in  the  source  region  Thus  the  result  of  the  depth  integration 
is  an  infinitesimal  dipole  which  has  the  form 

R(*J)  =  (2(u>/c)cos(0)Z?V))2  .  (16a) 

where  ^  ^ 

D'l")  =  h,;/)(y")*dy")  (16b) 
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It  follows  that  for  low  frequencies  the  cross-spectrum  p(u>;  z,  z',y)  is  approximately 
that  of  the  Cron  and  Sherman  model  for  m  =  1.  For  high  frequencies  the  dipole 
radiation  pattern  for  a  bubble  at  a  particular  depth  has  multiple  lobes  that  occur 
in  directions  which  vary  with  the  depth  of  the  bubble.  Thus  the  dipole  radiation 
pattern  for  a  bubble  at  one  depth  has  lobes  that  occur  in  the  directions  of  the 
nulls  in  the  radiation  pattern  for  a  bubble  at  another  depth.  The  effect  of  the  depth 
integration  in  (I4b)  is  to  produce  a  radiation  pattern  that  is  approximately  constant 
over  all  angles  for  which  the  multiple  lobes  occur.  This  can  be  seen  in  more  formal 
terms  by  rewriting  (14b)  in  the  form 


(u>;  y")  cos  (2(u>/ c)  cos(0)y")  dy" 


(17) 


For  high  frequencies  the  second  term  in  (17)  will  be  negligible  provided  that 
S,  b  (w;y")  varies  slowly  with  depth  in  comparison  to  cos  (2(u>/c)  cos(0)y").  In  the 
case  where  5tjj,(u>;  y")  »*  independent  of  y",  the  second  term  will  be  negligible  for  fre¬ 
quencies  and  angles  satisfying  /  >  (c/2Z))/  cos(0')  and  0  <  01 .  For  these  frequencies 
and  angles,  R( u>,  0)  ss  2.  Thus,  for  sufficiently  high  frequencies,  the  radiation  pattern 
is  approximately  that  of  a  .requency- independent  monopole  (except  for  angles  near 
£ir),  and  hence  the  cross-spectrum  p(u\  z,  z1,  7)  for  a  =  0  is  approximately  that  of  the 
Cron  and  Sherman  model  for  m  =  0.  For  later  reference  we  note  that  the  integrated 
spectrum  Cj(u>;  2,0,7)  for  a  =  0  has  the  frequency  dependence  (2(u >/c)D'(w))2  at 
low  frequencies  and  is  essentially  independent  of  frequency  at  high  frequencies. 


It  is  important  to  emphasise  that  the  simplified  models  are  obtained  with  the  source 
region  located  below  a  flat  surface.  The  bubbles  that  produce  the  sound,  however, 
are  assumed  to  be  located  well  within  the  breaking  wave  itself,  and  over  the  fre¬ 
quency  band  considered,  the  wavelength  of  the  radiated  sound  from  each  bubble  is 
comparable  to  or  smaller  than  the  dimensions  of  the  breaking  wave.  Consequently 
the  cross-spectra  obtained  from  the  simplified  models  might  differ  significantly  from 
those  obtained  if  the  shape  of  the  breaking  wave  were  taken  into  account. 
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4.  The  excitation  postulates 

In  Sect.  2  it  was  noted  that  for  linear  bubble  vibrations  excited  at  formation,  high- 
frequency  sound  is  emitted  only  if  the  bubbles  are  formed  in  a  non-equilibrium  state 
or  if  the  bubbles  are  subjected  to  a  rapid  change  in  pressure.  As  stated  in  Sect.  2, 
we  consider  only  those  bubbles  formed  as  the  whitecap  front  advances,  over-running 
a  layer  of  air.  In  this  section,  we  postulate  two  different  mechanisms  by  which  the 
bubble  vibrations  are  excited. 

The  first  postulate  is  essentially  a  pressure  excitation  of  the  bubbles.  We  assume 
that  at  formation,  the  bubbles  are  subjected  to  a  rapid  pressure  increase  arising  from 
the  weight  of  the  water  that  over-runs  the  bubbles.  The  over-running  is  assumed  to 
be  sufficiently  rapid  that  the  pressure  increase  approximates  to  a  step  increase.  The 
pressure  excitation  spectrum  is  then  given  by 

F(w;y,  R)  =  -ipgd(y)/u>,  (18a) 

and  the  excitation  energy  spectrum  is 

5f(u>;  y;  R)  =  (pgd(y)/u)\  (18b) 

where  d(y)  is  the  depth  of  the  bubble  formed  at  the  position  y. 

The  second  postulate  is  essentially  a  volume-rate  excitation.  At  the  moment  a  bubble 
is  formed  by  the  closure  of  an  air  cavity  we  assume  that  the  inward  velocity  of  the 
water  surrounding  the  cavity  imparts  an  initial  volume  rate  V'(O)  to  the  bubble. 
To  determine  the  initial  volume  rate  we  use  the  fact  that  for  the  small  bubble 
vibrations  for  which  the  linear  theory  is  valid  V(0)  %  4nR7R(0),  where  fl(0)  is  the 
initial  radial  velocity.  By  continuity,  R{0)  is  equal  to  the  inward  velocity  of  the 
water  at  the  moment  of  closure,  U.  Thus  by  identifying  jS(0)  with  V  we  obtain 

F(ur,y,R)  = -pRU.  (19a) 

To  obtain  the  excitation  energy  spectrum  we  assume  that  the  inward  velocity  at 
closure,  17,  is  a  random  variable  with  a  mean  that  is  independent  of  the  bubble 
radius  and  that  for  any  two  bubbles,  the  inward  velocities  at  closure  are  mutually 
uncorrelated.  These  assumptions  are  equivalent  to  assuming  that  at  the  moment  of 
closure  of  a  cavity,  the  velocity  of  the  surrounding  water  is  independent  of  the  size 
of  the  cavity  and  that  the  bubbles  are  formed  with  differences  in  position  and  time 
that  are  large  compared  with  the  correlation  distance  and  the  correlation  time  of 
the  turbulent  motion  in  the  front.  With  these  assumptions  the  excitation  energy 
spectrum  is  given  by 

Sf(w;  y,R)--=(pRV')\  (19b) 

where  U'  is  the  rms  inward  velocity  at  closure. 
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5.  Parameter  estimation  and  model  results 

In  this  section  we  present  examples  of  the  mean  bubble  spectrum  5m,b.(w)  and  the 
noise  intensity  spectrum  S(u>;z)  =  C(io\x,x),  for  the  two  excitation  postulates  de¬ 
scribed  in  the  previous  section.  The  examples  were  obtained  using  the  simplified 
model  for  the  ideal  medium  with  the  additional  assumption  that  the  bubble  gen¬ 
eration  rate  p(y),  the  pressure  transfer  spectrum  of  a  bubble  H(ur,y,  R),  and  the 
bubble-radius  density  p(R  |  y)  are  independent  of  the  bubble  position  y  throughout 
the  source  region.  The  simplified  model  equations  for  the  two  excitation  postulates 
are  derived  with  this  additional  assumption  in  Appendix  E. 

The  model  parameters  used  for  the  examples  were  estimated  for  a  wind  speed  of 
10  ms~*.  The  estimate  of  the  product  A  Mb  and  the  form  of  the  bubble-radius 
density  p(R)  were  based  on  a  model  for  the  generation  and  distribution  of  bubbles 
at  the  ocean  surface  proposed  by  Crowther  [15|.  In  this  model  the  bubbles  generated 
at  the  ocean  surface  are  distributed  by  the  action  of  turbulent  motions  that  carry 
them  downwards  and  by  buoyancy  that  forces  them  to  rise.  Under  steady-state 
conditions  the  number  of  bubbles  per  unit  volume  of  radius  R  that  are  distributed 
immediately  below  the  surface  is  given  by 

n(R)  =  s(R)/v(R),  (20) 

where  s(fi)  is  the  number  of  bubbles  per  unit  area  per  unit  time  of  radius  R  that  are 
generated  at  the  surface  and  v(R)  is  the  rise  velocity.  From  a  dimensional  argument, 
Crowther  concludes  that  the  surface  generation  rate  s(/()  is  proportional  to  (Wj R)3, 
where  W  is  the  wind  speed.  Using  this  relationship  in  (20)  with  an  expression  for 
v(f£),  and  fitting  the  resulting  bubble-radius  distribution  n(i£)  to  backscatter  data, 
he  concludes  that 

a(JZ)  =s  2.3  x  10-loW'r3f?-3.  (21) 

Using  this  model  an  expression  for  AMb  was  obtained  as  follows.  Since  A  specifies 
the  mean  number  of  breaking  waves  that  occur  per  unit  area  of  the  surface  per  unit 
time  and  Mb  specifies  the  mean  number  of  bubbles  generated  in  a  breaking  wave, 
AMb  specifies  the  mean  number  of  bubbles  generated  per  unit  area  of  the  surface  per 
unit  time.  Assuming  that  the  only  significant  source  of  surface-generated  bubbles  is 
breaking  waves,  AAfb  is  given  by  the  integral  of  s{R)  with  respect  to  f?,  i.e. 

t 

AMb  a  2.3  x  10",oVU3  /  R~3  d/J,  (22) 

where  Am|n  and  are  the  minimum  and  maximum  radii  of  the  bubbles  generated 
in  the  breaking  waves.  According  to  (22)  the  mean  number  of  bubbles  generated  in 
a  breaking  wave  is  proportional  to  the  integral  of  R~3  with  respect  to  R  and  hence 
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the  fraction  of  these  bubbles  with  radii  in  dR  is  proportional  to  R  3 .  It  follows  that 
the  bubble-radius  density  is  also  proportional  to  R~3 ,  i.e. 

/Amu 

p(R)*R-3/  R-3dR  (23) 


The  values  of  i?mln  and  Rm*x  determine  the  band  of  natural  bubble  frequencies  in  the 
mean  bubble  spectrum.  The  value  of  RmKX  was  inferred  from  laboratory  measure¬ 
ments  of  the  distribution  of  bubble  radii  in  simulated  breaking  waves  by  Monahan 
and  Zietlow  [16].  Over  the  range  of  radii  from  0.6  mm  to  5  nun,  their  distribution 
shows  an  R~3  dependence  -  the  same  as  Crowther’s  distribution  n(R)  for  a  rise 
velocity  that  is  independent  of  bubble  radius  (see  (20)  and  (21)).  Some  bubbles 
with  radii  between  5  and  9  nun  were  also  observed  in  the  laboratory  measurements. 
However,  the  possibility  cannot  be  excluded  that  such  large  bubbles  resulted  from 
the  coalescence  of  smaller  bubbles  formed  at  some  earlier  stage.  In  view  of  this  pos¬ 
sibility  we  restricted  ff.max  1°  5  mm,  corresponding  to  a  natural  bubble  frequency 
of  0.66  kHz.  For  Rm |n  we  used  a  value  of  0.3  mm  corresponding  to  10  kHz.  (We 
note  that  the  values  of  Rmin  and  Rmtx  also  affect  the  mean  bubble  spectrum  level 
owing  to  the  normalization  in  (23).  The  noise  intensity  spectrum  level,  however,  is 
not  affected  by  these  values  since  the  integral  in  (22)  and  its  inverse  in  (23)  cancel 
out.) 

An  estimate  of  A  was  obtained  from  A  %  w/AmTm,  where  to  is  the  fraction  of 
whitecap  coverage,  and  Am  and  Tm  are  the  mean  area  and  mean  lifetime  of  an 
individual  whitecap  event.  The  fraction  of  whitecap  coverage  was  taken  from  the 
observations  of  Ross  and  Cardone  [17].  In  their  definition  of  whitecap  coverage, 
these  authors  include  ‘large  new  foam  patches’.  This  was  taken  into  account  in 
estimating  the  mean  area  and  lifetime  of  an  individual  event.  Using  the  observation 
w  =  0.01  for  a  wind  speed  of  10  m  s-1 ,  with  estimates  Am  —  10  m*  and  Tm  =  10  s, 
we  obtained  A  =  10-4  breaking  waves  per  square  metre  of  ocean  surface  per  second. 

Using  the  above  values  of  Rmjn,  Rm»x  and  A  in  (22)  yielded  Mb  =  1.27  x  104  bubble 
generations.  (We  note  that  for  a  rectangular  bubble  generation  region  0.2  m  (deep)  x 
0.2  mx  2  m,  moving  at  constant  rate  during  a  breaking-  wave  lifetime  of  3  s,  1.27 x  104 
bubble  generations  implies  a  bubble  generation  rate  /»  -  5  3  x  104  bubbles  per  cubic 
metre  per  second.  Thus  in  a  cubic  decimetre  of  this  moving  region  53  bubbles  are 
generated  during  each  second  of  the  breaking- wave  lifetime.) 

Finally  for  both  excitation  postulates,  we  used  a  source  region  depth  D  =  0.2  m, 
and  for  the  volume-rate  postulate,  we  used  an  rms  inward  velocity  V  =  7  cm  s-1. 

The  mean  bubble  spectra  for  the  two  excitation  postulates  are  illustrated  in  Fig. 
3.  Both  spectra  have  their  highest  energy  levels  in  the  frequency  band  0.66-10 
kHz.  This  is  the  band  of  natural  bubble  frequencies  corresponding  to  the  range  of 
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Frequency  (Hz) 

Fig.  3.  Mean  bubble  spectra  for  0.2  m  source  region  depth:  dashed  curve 
-  pressure-jump  excitation;  dotted  curve  -  volume-rate  excitation  for  an 
rms  inward  velocity  of  7  cm  s“ 1 . 


bubble  radii  used.  Within  this  band  the  spectra  have  comparable  levels  at  the  lower 
frequencies  (large  bubbles)  and  both  decrease  in  level  as  the  frequency  increases 
(the  bubble  radius  decreases).  The  'pectra  differ  from  one  another  in  their  high- 
frequency  slopes.  For  the  pressure-jump  excitation  the  slope  is  about  -6  dB/octave 
over  the  full  band,  whereas  for  the  volume-rate  excitation  the  slope  varies  from 
about  -6  dB/octave  at  the  lower  frequencies  to  about  -2  dB/octave  e.t  the  higher 
frequencies. 

To  explain  the  slopes  of  the  mean  bubble  spectra  we  first  compare  the  corresponding 
single-bubble  spectra  5b(w;y,ii).  The  shapes  of  the  single-bubble  spectra  can  be 
seen  through  the  plots  of  Fig.  4.  These  plots  show  the  energy  transfer  spectrum 
|FT  (w;  -ft)!1  and  the  excitation  energy  spectra  Sjr(u/;  y ,  R )  for  a  3.3  mm  radius  bubble 
(/o  =  1  kHz),  at  the  rms  source  depth;  the  parameter  values  are  the  same  as  those 
for  the  mean  bubble  spectra  of  Fig.  3.  The  shapes  of  the  single-bubble  spectra 
are  seen  by  adding  the  corresponding  excitation  spectrum  to  the  transfer  spectrum. 
For  both  excitation  postulates  this  results  in  a  single-bubble  spectrum  with  a  sharp 
peak  at  the  natural  bubble  frequency  /o .  For  the  pressure- jump  excitation  the  single- 
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Fig.  4.  Single-bubble  spectrum  components  for  R  —  3.3 
mm:  solid  curve  -  energy  transfer  spectrum  at  1  m;  dashed 
line  excitation  energy  spectrum  for  pressure-jump  excita¬ 
tion;  dotted  line  -  excitation  energy  spectrum  for  volume-rate 
excitation. 


bubble  spectrum  approaches  zero  for  frequencies  on  either  side  of  /0  due  to  the  -6 
dB/octave  slope  in  the  excitation  spectrum.  In  contrast  the  single- bubble  spectrum 
for  the  volume-rate  excitation  only  approaches  zero  for  frequencies  less  than  /0  since 
the  excitation  spectrum  is  flat  -  for  frequencies  greater  than  /o  the  single-bubble 
spectrum  approaches  a  level  that  is  proportional  to  R*.  (We  note  that  the  peak 
source  level  of  the  single-bubble  spectrum  for  both  excitation  postulates  is  roughly 
90  dB  re  i  pPa  at  1  m.  This  level  results  from  an  amplitude  of  vibration  of  less  than 
1%  of  the  equilibrium  radius  which  is  well  within  the  linear  range.) 

Having  noted  the  difference  in  the  shapes  of  the  single-bubble  spectra  we  point 
out  that  for  both  excitation  postulates  the  peak  value  of  the  single-bubble  spec¬ 
trum  is  proportional  to  R*6~7.  This  can  be  seen  as  follows,  recalling  that  the 
transfer  spectrum  |/f(w;  ii)|J  has  a  peak  of  (R/6)*.  Por  the  volume-rate  excitation, 
Sf(w\ y,  R)  a  R1  so  that  the  peak  value  of  the  single-bubble  spectrum  is  propor¬ 
tional  to  For  the  pressure-jump  excitation,  Sf(ur,y,  R)  <x  u.--7  so  that  the 

peak  value,  which  is  well  approximated  by  |£T(woi  fi)|,5jr(wo;y,  A),  is  proportional 
to  (R/6 )*wc~2  or  to  RA6~7 .  For  natural  bubble  frequencies  less  than  10  kHz  (R  >  0.3 


% 


-  17  - 


SACLANTCEN  SR- 145 


mm)  S  is  only  weakly  dependent  on  R  (see  Eqs.  (3)),  and  thus  the  peak  value  for 
both  excitations  is  approximately  proportional  to  R 4 


Having  compared  the  single-bubble  spectra  the  slopes  of  the  mean  bubble  spectra  are 
explained  as  follows.  For  the  pressure-jump  excitation  the  single-bubble  spectrum 
goes  to  zero  for  frequencies  on  either  side  of  the  sharp  peak  at  u>o  so  that  for  a  fixed 
frequency  u>  the  major  contribution  to  the  integral  of  (8c)  occurs  when  u>o  is  approx¬ 
imately  equal  to  u.  This  occurs  for  bubble  radii  R  ^  2jr(3.3m  s-1)/u>.  For  these 
radii  the  integrand  is  approximately  equal  to  the  peak  value  of  the  single-bubble 
spectrum  weighted  by  the  bubble-radius  density.  The  peak  value  is  proportional  to 
R*  and  the  radius  density  is  proportional  to  R~3.  Consequently  the  integrand  is 
proportional  to  R  so  that  the  integral  is  proportional  to  R 2  or  equivalently  to  u>~ 2 
in  accordance  with  the  -6  dB/octave  slope  seen  in  the  mean  bubble  spectrum  of 
Fig.  3. 

For  the  volume- rate  excitation  the  preceding  argument  does  not  apply,  even  though 
the  peak  of  the  single-bubble  spectrum  has  the  same  radius  dependence  as  that 
for  the  pressure-jump  excitation.  This  is  because  the  single- bubble  spectrum  does 
not  go  to  zero  for  frequencies  greater  than  a»o,  rather  it  approaches  a  level  that  is 
proportional  to  R* .  As  a  result,  when  w  is  large  contributions  to  the  integral  of  (8c) 
from  the  single-bubble  spectra  for  small  u>o  (large  R)  are  not  small  in  comparison 
to  those  for  u>o  «  u>  and  hence  cannot  be  neglected.  In  tact,  it  is  the  integration 
of  these  contributions  that  results  in  the  changing  slope  at  high  frequencies  in  the 
mean  bubble  spectrum  of  Fig.  3. 

The  examples  of  the  noise  intensity  spectra  for  the  two  excitation  postulates  are 
shown  in  Fig.  5.  Also  shown  is  the  Wenz  spectrum  for  the  10  m  s~’  wind  speed 
taken  from  (1).  Both  model  spectra  show  a  near-constant  level  at  the  lower  fre¬ 
quencies  of  the  band  0.66-10  kHz,  followed  by  a  decrease  in  level  at  the  higher 
frequencies.  At  the  lower  frequencies  the  near-constant  level  in  both  spectra  results 
from  the  cancellation  of  the  dependence  in  the  mean  bubble  spectrum  by  the 
dependence  in  the  integrated  propagation  spectrum  (see  Sect.  3).  At  the  higher 
frequencies  the  slopes  of  the  two  model  spectra  are  approximately  those  of  the  cor¬ 
responding  mean  bubble  spectra.  This  results  from  the  fact  that  at  high  frequencies 
the  integrated  propagation  spectrum  is  essentially  independent  of  frequency.  (For 
frequencies  greater  than  about  7  kHz,  the  integrated  propagation  spectrum  decreases 
with  frequency  because  of  the  absorption  loss  -  computed  for  a  hydrophone  depth 
of  1000  m  -  and  hence  the  model  noise  intensity  spectra  fall  off  more  rapidly  than 
the  corresponding  mean  bubble  spectra. ) 


The  noise  intensity  spectrum  for  the  pressure-jump  excitation  shows  better  agree¬ 
ment  with  the  Wenz  spectrum  than  that  for  the  volume-rate  excitation  at  the  high 
frequencies.  At  these  frequencies  the  slope  of  the  pressure-jump  spectrum  is  only 
.slightly  more  negative  than  the  -5  dB/octave  slope  of  the  Wenz  spectrum,  whereas 
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Fig.  5.  Noise  intensity  spectra  for  10  in  s"1  wind  speed:  solid  <  ^rve 
-  Wenz  spectrum;  dashed  curve  -  pressure-jump  excitation  and  dotted 
curve  -  volume-rate  excitation  for  A  =  10  ~ 4  breaking  waves  m"1  s_1, 
p(R)  oc  R~3,  R  between  0.3  and  5  min,  —  1.27  x  10*  bubbles. 

the  slope  of  the  volume-rate  spectrum  is  at  least  2  dB  more  positive  than  that  of 
the  Wenz  spectrum. 
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6.  Summary  and  conclusion 

A  general  model  has  been  developed  for  the  generation  of  high-frequency  noise  by 
linear  vibrations  of  the  bubbles  formed  in  breaking  waves.  This  model  essentially 
comprises  two  cross-spectra:  the  breaking-wave  cross- spectrum  that  describes  the 
noise  produced  by  individual  breaking  waves,  and  the  noise  cross-spectrum  that 
describes  the  noise  produced  by  the  aggregate  of  the  breaking  waves.  The  breaking- 
wave  cross-spectrum  is  obtained  by  assuming  that  the  noise  from  the  aggregate 
of  bubbles  is  the  superposition  of  the  noise  from  individual  bubbles  and  that  the 
bubbles  are  generated  according  to  a  Poisson  process.  This  cross-spectrum  can  be 
written  as  the  sum  of  two  cross-spectra,  one  of  which  makes  the  dominant  energy 
contribution  at  high  frequencies.  The  noise  cross-spectrum  is  obtained  by  assum¬ 
ing  that  the  pressure  waveforms  from  different  breaking  waves  are  statistically  in¬ 
dependent  and  identically  distributed,  and  that  the  breaking  waves  are  generated 
according  to  a  second  Poisson  process. 

Two  simplified  models  have  been  obtained  from  the  general  model  by  assuming  that 
the  source  region  lies  below  a  flat  surface  and  that  the  propagation  is  determined 
as  that  for  either  a  stratified  medium  or  an  ideal  medium.  As  such,  these  models 
can  be  viewed  as  extensions  of  the  Kuperman  and  Ingenito  model  and  the  Cron  and 
Sherman  model  for  the  case  of  noise  generated  by  bubbles  in  breaking  waves.  Both 
simplified  models  characterize  the  noise  source  in  terms  of  a  spectrum  at  each  point 
on  the  ocean  surface.  For  the  stratified  medium  the  noise  cross-spectrum  is  obtained 
from  a  depth  integral  of  the  ‘source-sheet’  cross-spectrum  of  the  Kuperman  and  In¬ 
genito  model.  For  the  ideal  medium  the  noise  cross- spectrum  is  obtained  using  a 
frequency-dependent  radiation  pattern  in  the  Cron  and  Sherman  model. 

Examples  of  the  noise  intensity  spectra  for  a  10  m  s'1  wind  speed  have  been  pre¬ 
sented  for  two  postulated  excitation  mechanisms  -  a  pressure-jump  excitation  and 
a  volume-rate  excitation.  These  spectra  were  computed  using  the  simplified  model 
for  the  ideal  medium  with  values  of  the  parameters  mostly  inferred  from  published 
measurements.  Although  different  the  two  spectra  both  show  levels  comparable  to 
the  Wenz  spectrum  for  the  same  wind  speed.  Of  the  two  spectra,  the  high-frequency 
slope  of  the  pressure-jump  spectrum  is  in  better  agreement  with  the  Wenz  spectrum. 
It  is  noted,  however,  that  the  idealized  propagation  assumption  used  in  the  simpli 
fied  model  neglects  the  shape  of  the  breaking  wave.  If  this  shape  were  accounted  for 
there  would  be  a  different  frequency  dependence  in  the  radiation  pattern  which  could 
significantly  change  the  shape  of  the  model  spectra  at  the  higher  frequencies.  It  is 
also  noted  that  there  are  uncertainties  in  the  overall  levels  of  the  model  spectra  due 
to  uncertainties  in  the  parameters,  Nevertheless,  the  rough  agreement  between  the 
model  spectra  and  the  Wenz  spectrum  supports  the  hypothesis  that  linear  vibrations 
of  bubbles  generated  in  breaking  waves  are  a  primary  source  of  wind-related  noise. 
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Appendix  A 

The  breaking-wave  and  bubble  occurrences 


In  this  appendix  we  obtain  the  multi- dimensional,  non- homogeneous  Poisson  pro¬ 
cesses  for  the  breaking-wave  occurrences  and  the  bubble  occurrences.  To  this  end 
we  first  summarize  the  relevant  properties  of  a  general  Poisson  process. 

Points  are  said  to  be  distributed  in  a  space  „Y  according  to  a  non-homogeneous 
Poisson  process  if  for  any  subset  A  contained  in  X  the  number  of  points  in  A  is  a 
Poisson  random  variable  N  with  the  probability  that  n  points  lie  in  A  given  by 

Pt[N  =  »)  =  e-MMn/n\,  (A.l) 

where 

M  —  J  A(x)  d3x  (A. 2) 

is  the  mean  number  of  points  in  A  and  A(x)  is  the  rate  function  (mean  number 
of  points  per  unit  volume)  of  the  process.  Clearly,  to  specify  a  non-homogeneous 
Poisson  process  in  the  space  X  it  suffices  to  specify  its  rate  function  A(x).  In  the 
special  case  where  the  rate  function  is  independent  of  x,  the  Poisson  process  is  said 
to  be  homogeneous.  In  this  case,  M  is  given  by 


M  =  Avolume(A).  (A.3) 

Given  that  there  are  n  points  in  A ,  the  Poisson  process  describes  those  points  in 
terms  of  a  joint  probability  density  given  by 


p(x1(...,x„)  =  J|p(xfc), 

fc= 1 


(A. 4a) 


where 


A(x)/M, 

0, 


x  6  4 
otherwise. 


(A.4b) 


It  follows  from  Eqs.  (A. 4)  that  the  points  are  independent,  identically- distributed 
random  variables  with  the  common  probability  density  of  (A. 4b). 


The  Poisson  process  for  the  breaking- wave  occurrences  is  obtained  by  identifying 
the  point  x*  with  (rje.r*)  and  taking  A  to  be  of  the  form  A  =  A%  X  ©,  where  A, 
is  the  total  area  of  the  ocean  surface  where  breaking  waves  can  occur  and  0  is  the 
time  interval  [-|T,  \T\.  It  is  assumed  that  during  the  time  interval  of  interest  the 
rate  function  depends  only  on  the  breaking-wave  occurrence  position  r  and  not  on 
the  occurrence  time.  The  number  of  breaking  waves  generated  in  the  surface  region 
A,  during  the  time  interval  0,  Nw,  has  the  probability  function  of  (A.l)  with  mean 


Afw 


(A. 5) 
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where  the  units  of  A(r)  are  whitecap  occurrences  per  unit  area  per  unit  time.  Fur¬ 
thermore,  given  that  there  are  n  breaking  waves  in  A,,  the  occurrence  positions  and 
times  of  those  breaking  waves  are  described  in  terms  of  the  joint  probability  density 


where 


p(*T,n, 

•  ■ • i rn>  Tn) 

*=1 

(A.6a) 

P(r,r) =  | 

A(r)/Mw, 

.0, 

(r,r)  €  A,  x  © 
otherwise. 

(A.6b) 

The  Poisson  process  for  the  bubble  occurrences  is  obtained  by  identifying  x*  with 
(yfc,  Vk)  and  taking  A  to  be  of  the  form  A  —  V  x  6y,  where  the  bubble  source  region 
V  is  the  total  region  swept  out  by  the  bubble  generation  region  during  the  lifetime 
of  the  breaking  wave  and  Qy  is  the  time  interval  during  which  the  point  y  £  V  is 
contained  in  the  generation  region.  Thus,  A  represents  the  space-time  set  in  which 
the  bubbles  are  generated  and  hence  the  rate  function  n(y,v)  is  non-zero  only  for 
(y,  v)  €  A.  The  number  of  bubbles  generated  in  the  source  region  during  the  lifetime 
of  the  breaking  wave,  iVb,  has  the  probability  function  of  (A.l)  with  mean 

Mb  =  /  /  n(y<v)  di/d3y,  (A.7) 

Jv 

where  /i(y,  v)  has  the  units  of  bubble  occurrences  per  unit  volume  per  unit  time. 
Given  that  there  are  n  bubbles  in  A,  the  occurrence  positions  and  times  of  those 
bubbles  are  described  in  terms  of  the  joint  probability  density 


p(yi  >  v\ ,  •  •  • ,  y«,  vn)  =  p(y*,  v k ), 

1 


(A. 8a) 


where 


(y,  v)  €  V  x  ©y 
otherwise. 


(A. 8b) 


Using  (7a),  the  single-bubble  occurrence  density  can  also  be  written  in  the  form 


PiY,v)  =  Afb  V(y 1  y),  (A.9a) 

where 

I  y )  =  M(y.  *0/#*(  y)  (A.9b) 

is  the  conditional  occurrence-time  probability  density.  This  probability  density  de¬ 
scribes  the  bubble  occurrence  times  at  a  fixed  point  within  the  total  volume  V. 

Finally,  from  (A. 9a)  and  the  assumption  that  the  bubble  radius  depends  only  on 
the  occurrence  position  of  the  bubble  and  not  on  its  occurrence  time,  the  joint 
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probability  density  of  the  bubble  occurrence  position,  its  occurrence  time  and  its 
radius  can  be  written  as 

n 

p(yi,  yn,  *%,  Rn)  =  fj  p{yt.,Vk,Ri,),  (A.iOa) 

fc=i 

where 

p(y,^,  R)  =  Afb”  l/‘(y)/‘(«/ 1  y)p(R  I  y)-  (A.iob) 
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Appendix  B 

Derivation  of  the  noise  cross-spectrum 

Let  p(t;z)  be  the  total  pressure  observed  on  a  hydrophone  located  at  position  z 
due  to  the  noise  generated  by  the  aggregate  of  breaking  waves.  We  assume  that  the 
contributions  from  the  individual  breaking  waves  are  statistically-iudependent,  zero- 
mean  processes  with  a  common  probability  law  that  is  independent  of  the  position 
and  the  time  at  which  the  breaking  wave  occurs.  Under  this  assumption,  p(t;x)  can 
be  written  in  the  form 

p(t\z)  =  ^Tp*(t  -  rfc;rfc,z),  (B.l) 

where  pw(t  -  rfc;rk,z)  is  the  pressure  contribution  from  the  Ifeth  breaking  wave  to 
occur,  rfc  is  the  time  at  which  the  breaking  wave  begins,  and  r*  is  its  position  on  the 
ocean  surface.  The  noise  cross-spectrum  between  the  pressures  observed  at  positions 
zi  and  z2  is  defined  as  the  fourier  transform 

C(u>;z,,z2)  =  /  R(t,t  +  u;zi,zi)e  ,wudu,  (B.2a) 

J  -oo 

where 

R{t,t  +  u;xt,z2)  ^  E[p(f;  z,  )p(t  +  u;z2)]  (B.2b) 

is  the  space-time  correlation  function  between  the  pressures  at  positions  zt  and  Z2 
and  E[-}  is  the  expected  value  operator  over  all  stochastic  variables  determining  the 
received  pressure.  Note  that  the  noise  cross-spcctrum  is  well  defined  only  when  the 
space-time  correlation  function  does  not  depend  on  the  time  t. 

In  this  appendix  we  show  that  if  the  breaking- wave  occurrence  positions  and  oc¬ 
curence  times  {(r*.,Tfc)}  are  described  by  the  Poisson  process  specified  in  Eqs.  (A. 5) 
and  (A.6)  then  the  noise  cross-spectrum  is  given  by 


C(w; xi , *j)  =  /  A(r)(?w(w; r, r.\ , z2) d*r,  (B.3a) 

Ja. 

where 

C'w(w;r,zi,zJ)  =  £w[p;(W;r,zt)Pw(w;r,z2)]  (B.3b) 

is  the  breaking-wave  cross-spectrum,  £w[]  is  the  expected  value  operator  for  fixed 
breaking-wave  occurrence  positions  and  occurrence  times,  and 


Pw(u>;  r,Zj) 


/; 


pw(t;r,Zj)e 


-iu>t 


dt 


(B.3c) 


is  the  pressure  spectrum  resulting  from  an  individual  breaking  wave.  The  noise 
intensity  spectrum  is  obtained  by  setting  Z)  -  z2  =  z  in  Eqs.  (B.3). 


» 
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To  establish  Eqs.  (B.3)  it  suffices  to  show  that 


where 


+  =  /  A(r)JZw(ti;r,*i,«i)d*r,  (B.4a) 

Ja. 

JJw(u;r,*i,«a)  =  /  £w[pw(t;r,*i)pw(<  + u;r,*i)]dt.  (B.4b) 

J  00 


Given  this  result,  (B.4a)  can  be  used  in  (B.2a)  to  write  the  noise  cross-spectrum  in 
the  form  of  (B.3a)  where 


Cw(u> 


=  J 


flw(u;r,*i,*j)e  ,u,udu. 


Eqs.  (B.3b)  and  (B.3c)  follow  from  using  (B.4b)  in  (B.5). 

To  derive  Eqs.  (B.4),  we  first  obtain  an  expression  for  the  joint  characteristic  function 
of  the  random  variables  p(t;  >1 )  and  p(t  -f  u;  a2 ) 

*(«,,«!)  =  j B[c‘(''tP(‘i-i)+vji<‘+«i»*))]i  (B.6) 

and  then  use  the  moment  theorem 

£[p(t;zl)p(t  +  u;*i)j  =  -d*$(0,0)/dDi3u,  (B.7) 

to  determine  the  expected  value  in  (B.2b).  For  notational  convenience,  we  write 
p(t\  *1)  and  p(t  +  u;  *j )  in  the  form 


where 


p(t;*i)  =  ^Pfc,i(rfc,rfc), 
k 

p(t+  u;»2)  =  ^Pfc.a(r),,  r»), 


Pfc.l (*"fc .  Tfc)  =  Pw(<  -  Tfc \rk,*i), 

I,rfc)  =  pw(t  +  u  -  Tfc;  rj,,  «j). 


(B.8a) 


(B.8b) 


(B.9a) 

(B.9b) 


The  characteristic  function  is  obtained  by  expanding  the  expected  value  in  (B.6)  in 
terms  of  expected  values  conditioned  on  the  event  that  there  are  exactly  n  breaking 
waves  in  the  region  A  =  A,  X  0.  The  result  is 

OO 

*(»,, v,)  =  £  £[c«vj j  Nm  n)Pr(ATw  =  „),  (B.lOa) 
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where 

£;[e'(«iP(*i»i)+,,iMt+ui*j))  |  -  n]  = 

Ert.n . . .  >+v’p<‘+“1,,))]]  ,  (B.lOb) 

Pr(jV„  =  n]  is  given  by  (A.l)  with  M  given  by  (A.5),  and  E,x ,n ,„[•]  is  the 

expected  value  operator  over  the  breaking-wave  occurrence  positions  and  occurrence 
times.  Using  Eqs.  (B.8),  the  expected  value  in  the  nth  term  in  the  sum  can  be  written 
as 

. »n*1+w.»|] 

-  . 

=  a„ . . .  [fi-in 

h=l 

=  f[  Etk,rk  [£w[e<(v‘p*.»+,'lpM>]] 

kz si 

=  (E'k, rk  [£w(c<<u‘p*.‘+”lP*'»)]])n  ,  (B.ll) 


where  we  have  used  the  Fact  that  the  points  {(t*,7"*)}  are  statistically  independent 
and  identically  distributed  to  replace  the  expectation  of  the  product  by  the  product 
of  the  expectations  in  the  third  equation  and  to  replace  the  product  by  the  power 
in  the  last  equation.  Next,  we  write  the  outer  expected  value  in  the  last  equation  as 
an  integral  using  (A.6b)  for  the  probability  density.  The  result  is 

Er„,rk  [JEw[c<<WlPfc  l+Va,>*a>l]  = 

/r/i  , 

/  M"1A(r)£;w[e<t,'lP*.i+u»p'-.»>]dJrdr.  (B.12) 

-r/j  J a. 

Using  (B.12)  in  the  last  equation  of  (B.ll)  and  then  using  the  result  in  Eqs.  (B.10) 
and  substituting  for  Pr[JV„  =  n]  using  (A.l)  yields 


#(*i,t>,)  =  e'M«  £  ^  j*  M;lA(r)£w[  c<<”iPM+*m.J>]d,rdr  j 


x  M£/n\. 


(B.13) 


The  desired  expression  for  the  characteristic  function  follows  by  writing  the  sum¬ 
mation  as  an  exponential  function  and  using  (A.5)  to  express  Afw  in  terms  of  A(r). 
The  result  is 

*(«,,«*)  =  exp  ^  j ^  A(r)  (£.[c«*'‘p^+p’p‘.’>]  -  l)  d’rdrj  .  (B.14) 
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Eqs.  (B.4)  are  obtained  by  first  using  (E.14)  to  substitute  for  the  characteristic 
function  in  (B.7),  differentiating,  setting  u,  =b  vj  =  0  and  using  Eqs.  (B.9)  to 
substitute  for  p/,,i  and  p^j.  The  result  is 

E(p(i;ai)p(<  +  u;a,}j  = 

tT,i  r 

/  /  A(r)I?w[pw((  -  r;p,*i)pw(t  -(-  u  —  r;r,  Bj)]d,rdr.  (B.15) 

J-T/i  J A, 

Interchanging  the  order  of  integration,  making  the  change  of  variable  r'  =  f  -  r,  and 
taking  the  limit  as  T  approaches  infinity  results  in  Eqs.  (B.4). 
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Appendix  C 

Derivation  of  the  breaking-wave  cross-spectrum 

In  the  previous  appendix  we  derived  the  expression  for  the  noise  cross-spectrum  in 
terms  of  the  breaking- wave  cross-spectrum.  In  this  appendix  we  derive  the  expression 
for  the  breaking-wave  cross-spectrum  of  (B.3b)  under  the  assumptions  of  Sect.  2.  For 
generality  we  allow  the  pressure  transfer  spectrum  y,  R)  as  well  as  the  pressure 
excitation  spectrum  F(u;y,R)  to  be  stochastic  and  define  the  single-bubble  cross- 
spectruin  Cb(w;  y.y',  R,  R')  and  the  single-bubble  spectrum  Sb(u>;y,.R)  by 

CKwjy.y'.Jl,*')  =  Eb[P:(^y,R)Pb(^y\R')},  (C.la) 


and 

Sb(W;y,i*)=  £b(|Fi,(u,iy,fi)|2],  (C.lb) 

where  £b[  ]  is  the  expected  value  operator  for  fixed  bubble  occurrence  positions, 
occurrence  times  and  radii.  Note  that  when  the  pressure  transfer  spectrum  is  deter¬ 
ministic,  Cb(u»;y,y',  R,  R')  and  5b(u;;  y,  R)  are  determined  by  the  excitation  cross- 
spectrum  and  the  excitation  energy  spectrum  of  (6a)  and  (6b). 

We  show  that  under  the  assumptions  of  Sect.  2,  the  breaking-wave  cross-spectrum 
is  given  by 


C„(w;r,*i,**)  =  /  /i(y)h’r.b.(u>;y)(;*(w,r  +  y,*,)G(w,r  +  y,z2)d*y 
Jv 

*u  My)My')^r.b.(w;y,y')0‘(u>,r  +  y,*i)G(u>,r  +  y',z2) 

x  $*(w  |  y)$(w  |  y')d3y  d3y',  (C.2a) 

where 

S'r.b(w;y)=/  Sb(w,y,R)p{R  \  y)dR  (C.2b) 

Jo 

is  the  radius- averaged  bubble  spectrum, 

CVb.(u>;y,y')=  l°°  f  C’b(u»;y,y',  R,R')p(R  \  y)p(R'  \  y')dRdR'  (C.2c) 
Jo  Jo 

is  the  radius-averaged  bubble  cross-spectrum,  and 


*(w  I  y)  =  /  |  y)e,wl/dts 

Je, 


(C.2d) 


is  the  characteristic  function  of  the  bubble  occurrence-time  probability  density  ft(i/  | 

y). 


* 
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The  cross-spectrum  of  (C,2a)  is  derived  by  expanding  the  expected  value  operator 
in  (B.3b)  as  a  sum  of  conditional  expectations  and  then  evaluating  each  term  in  the 
sum.  The  expansion  has  the  form 


Cw(w;r, ■,,*»)  =  £Jw[P„(u>;r,zi)Pw(u>;r,*2)  I  Ab  =  njPr[^b  =  «].  (C.3a, 


n=l 


where 


|  Nh  =  n]  = 

[^b[Pw(w;  r,  Zi  )Pw(a;;  v,  *j)]] ,  (C.3b) 

Pr[ATb  =  n)  is  given  by  (A.l)  with  M  given  by  (A. 7),  and  Eyi>t,liHl . . ,*„[•]  is 

the  expected  value  operator  over  the  bubble  occurrence  positions,  occurrence  times 
and  radii.  To  simplify  the  notation,  we  rewrite  the  expression  for  the  spectrum  of 
the  received  pressure  (5)  in  the  form 


Nb 

-  £/Uy*.^)G>(y*)e,W\  (C.4) 

fc=i 

where  i\,(y,  R)  and  Gj( y)  are  the  bubble  pressure  spectrum  and  the  function  (?(u>, 
r-f  y ,  Zj)  for  fixed  values  of  the  frequency,  the  breaking- wave  occurrence  position  and 
the  hydrophone  positions.  Using  (C.4)  we  write  the  inner  expected  value  in  (0.3b) 
in  the  form 


£i,(p;(u-;r,zi)Pw(w;r,Z2))  - 


£  Eb[(Pb(yfc,^)<3l(yfc)c,wM'Pb(yfc.,Pfc0G2(yfc.)^').  (C.5) 

fc,fc'= i 


Rewriting  the  double  summation  in  (C.5)  as  a  sum  over  the  diagonal  terms  and 
a  sum  over  the  off-diagonal  terms  and  using  the  definitions  of  the  single-bubble 
cross-spectrum  and  the  single-bubble  spectrum  in  Eqs.  (C.la)  and  (C.lb)  results  in 


Eb[P,;(w;r,iii)Pw(u>;r,zj)]  =  5i>(w;yfc,  Rk)Gi(yk)G2(yk) 

fc=l 

n 

+  V]  Ch(w;yk,yk,,R  ,  Rk>)Gi(yk)Gi{yk,) 

k^k1 
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Using  (C.6)  in  (C.3b)  and  interchanging  the  expected  value  with  the  summation 
yields 


£w[/,;(u);r,z,)/>w(u;;r,Jij)  |  Nh  =  n] 

n 

=  E  ^n.n.Ri . yn.^n.Rr  [Sb(u>;yfc,  flfc)Gi*(yfc)G2(yfe)] 


fc=i 


+  E  . yn,K„.Rn[Cb(w;yfc,yfc..iZfc,flfc.)Gt*(yfc)G2(yfc.)c-^e<“1'*'] 

k*k' 


n 

-  E  £;yi..*'fc.Riil^b(u';  y>'  ^fc)Gi(y»«)G2(yfc)] 

k=i 

n 

+  E  •Ey«..‘'»i.Rii.yi1t.*'i,<.R|,<[Gb(^;yje,yfc',-Rfc,-Rfc')G?(yfc)G2(yfc>) 
k*h> 

x  e-w*e<w^'].  (C.7) 


Using  the  joint  probability  density  for  (y*,  t/k>  Rk)  of  (A. 10b),  the  first  expected 
value  in  (C.7)  can  be  written  as 


Eyfc.Kfc,R*[^b(w;yfc,fl)t)G1*(yfc)G!2(yfc)] 

=  j  A/b“V(y k)n{vk  I  yfe)p(R*  I  yfc)^b(w;yfc,  Rk)Gl(yk)G2(yk)d3yk  di/fc  dRk 

=  Mb-‘  /  p(y)5r.b,(w;y)<jJ(y)G2(y) d3y,  (C.8a) 

Jv 


where  5r,b,(w;y)  is  the  radius-averaged  bubble  spectrum  of  (C.2b).  Similarly,  the 
second  expected  value  in  (C.7)  can  be  written  as 


^yfclKfc.fifc,yM.»/*',^fc*(Gb(a;;yfc,yfc*,iifc,/Jfc.)<5,*(yfc)(72(yfc.)e  »"‘'fce,<*,‘V] 

=  // Afb'V(yfcM vk  I  y k)p{Rk  I  yfc)p(yfc')p(»/fc'  I  yfc< )p(-Rv  I  yk>) 


x  Cb(u;;  yk,  yk . ,  Rk,  Rk,  )G1*(yfc)G2(yfc. ) 

X  e  "•w‘'*e‘"‘V  d3yfc  d vk  d Rk  d3yfc»  dvk>  d-/2fc» 

=  p(y)My')Gr.b.(w;y,y')Gl*(y)G2(y') 

x  $*((1;  |  y)$(w  |  y')d3yd3y',  (C.8b) 


where  Cr  b^iytyO  i*  the  radius- averaged  bubble  cross-spectrum  of  (C.2c)  and 
$(w  |  y)  is  the  characteristic  function  of  (C.2d).  Using  (C.8a)  and  (C.8b)  in  (C.7) 
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and  evaluating  the  summation  yields 

£w[P;(w;r,z1)Pw(w;r,2i2)  |  Nb  =  n)  = 

nMbl  /  P(y)St.b.(w;y)Gi'(y)G2(y)d3y 

Jv 

+  n(n  -  l)Mb~J  f  (  Li(y)n(y')Ct.h.(w\y,y')Gul(y)Gi{y') 

Jv  Jv 

x  |  y)$(u>  |  y')d3y  d3y’.  (C.9) 

Finally,  using  (C.9)  in  (C.3a)  and  the  facts  that 

OO 

Mb  =  ^nPrflVb  =  n]  (C.lOa) 

n=1 

and  for  a  Poisson  random  variable 

OO 

Mb (Mb  +  1)  -  nJ  Pr[lVb  =  n)  (C.lOb) 

n=  1 

yields  (C.2a). 
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Appendix  D 

The  high-frequency  approximation 

The  breaking- wave  energy  spectrum  is  obtained  from  (C.2a)  by  setting  Z\  =  z2  =  z. 
The  result  can  be  written  in  the  form 


Sw(w;  r,  z)  =  5wi  (w;  r,  z)  +  Sw2(u>;  >*,  z), 


(D.la) 


where 


5wi(u>;r,z)  =  f  /t(y)5r.b.(w;  y)|(?(w, r  +  y,  z)|2  d3y,  (D.lb) 

Jv 

5w2(u>;r,z)  =  f  f  /u(y)^(y')Crb.(w; y , y')G*(w, r  +  y, z)G(u»,r  +  y', z) 

Jv  Jv 


X  $ '(w  I  y)$(u>  I  y')  d3y  d3y'. 


(D.lc) 


In  this  appendix  we  derive  a  frequency  dependent  upper  bound  on  the  ratio  of  the 
second  term  in  (D.la)  to  the  first  and  use  that  bound  to  identify  the  frequencies 
where  the  contributions  of  the  second  term  can  be  neglected. 

The  upper  bound  depends  only  on  the  mean  number  of  bubbles  generated  in  the 
breaking  wave  and  the  probability  density  of  the  bubble  occurrence  times.  Ar  such  it 
is  completely  determined  by  the  bubble  generation  rate  function  p(y,  u).  The  formal 
statement  of  the  upper  bound  is  as  follows.  Let  p(w)  be  defined  by 

p(w)  =  Afb|$(w  |  y0)|J,  (D.2a) 

where  d>(w  |  y)  is  the  characteristic  function  of  the  probability  density  n(u  |  y)  and 
y0(w)  is  determined  by 


|#(w  |  yo(w))|J  =  msuc|$(w  |  y)|2.  (D.2b) 

y  fV 

Then  the  spectra  5wi(w;r,z)  and  5w2(u';r,z)  satisfy 

|Sw2(w;r,z)|  <  p{w)6'wi(w;r,z).  (D.2c) 


To  apply  the  upper  bound,  we  consider  the  special  case  where  the  bubble  generation 
region  has  fixed  dimensions  and  moves  with  constant  speed.  We  further  assume  that 
the  bubbles  in  this  region  are  generated  at  a  constant  rate.  For  this  case,  n(y,v)  is 
constant  for  each  v  €  0y  and  @y  is  independent  of  y  e  V.  It  follows  from  (A. 9b) 
that  fi(i>  |  y)  is  a  uniform  probability  density  for  each  v  £•  0y .  Thus  by  using  the 
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expression  for  the  characteristic  function  of  a  uniform  random  variable  in  (D.2a)  we 
obtain 


Mb 

lh“>TY 


(D.3) 


where  T  is  the  duration  of  the  interval  @y.  For  Mb  —  10*  bubbles  and  T  =  0.1  s, 
p  <  0.1  for  frequencies  /  >  1  kHz.  Thus,  for  the  high  frequencies  of  interest,  the 
energy  spectrum  is  well  approximated  by  the  first  term  of  (D.la). 

The  derivation  of  the  upper  bound  makes  use  of  the  following  prop  arty  of  expected 
value  operators.  Let  A’  and  Y  be  two  complex  random  variables.  T  ->en 


jf?[AF]|  <  £{)AFj]  <  (JE[|AtT|2])1^2(JE7(|V|2])1/3  (D.4a) 


and 

E(|Ar|,/Jj  <  (JS[|JT|])1/a.  (D.4b) 

The  inequality  (D.4a)  is  simply  an  application  of  Cauchy’s  inequality  where  the 
inner  product  is  the  expectation  operator.  The  inequality  (D.4b)  is  obtained  from 
(D.4a)  by  setting  V  —  1  and  using  | Af  | 1  in  place  of  |X|. 

To  derive  the  upper  bound,  we  first  show  that  the  magnitude  of  the  radius- averaged 
cross-spectrum  satisfies 

|0't.b.(w;y,y')|  <  (5r.b.(w; y )Srb.(w;  y'))1  /2  .  (D.5) 


To  this  end  we  rewrite  (C.2c)  in  the  form 


C'r.b.(u>;y,y')  =  f  f  (Sb(u>; y,  fl)5b(w; y',  R'))l/7  Cb(w;y , y\ R,  R') 
J  0  Jo 


xp(R\y)p(R'\y')dRdR', 


(D.Sa) 


where 


Cb(w,y,y',R,R')  = 


Ch(u;yy,R,R') 

(Sh{w,y,R)Sh(w}y',R'))l/2 


(D.6b) 


is  the  normalized  single-bubble  cross-spectrum.  By  applying  (D.4a)  to  Cb(u>;y,y', 
R,  R' )  =  FJbffb^u’jy,  y',  i?' )j,  and  then  using  the  result  in  (D.6b),  we  obtain 


\Cb(w,y,y\R,R')\<  1. 


(D.7) 
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Taking  the  magnitude  of  both  sides  of  (D.6&)  and  using  (D.7)  results  in 

poo 

KVi>.(“';y,y')l  <  /  (5b(w;y,  fl))1/2p(fl  |  y)  dT. 

JO 

x  (°°  (5b(u»;  y\  R))i/2  p(R  |  y')  d R.  (D.8) 

J  o 

Applying  (D.4b)  to  both  integrals  on  the  rhs  of  (D.8)  with  |X|  equal  to  5b(w;y,  R) 
and  5b(w;y',  R)  respectively  we  obtain 

|Cr.b.(w;y,y,)|  <  (jH  Sh(u‘,y,R)p{R  \  y)dR)/2 

Of  CO  \  1  /2 

1  5b(w;y',i?)p(«|y')diiJ  .  (D.9) 

The  inequality  (D.5)  follows  immediately  from  (D.9)  and  the  definition  of  the  radius- 
averaged  bubble  spectrum  (C.2b). 

We  now  complete  the  derivation.  By  taking  the  magnitude  of  both  sides  of  (D.lc) 
and  using  (D.5)  we  obtain 

Sw2(a>;r,z)<  (  j  p(y)p(y')(‘S'r  b.(w;y)St.b.(^;y'))1/2  |G*(w,r  +  y,z) 

Jv  Jv 

x  G(u>,r  +  y',z)$*(w  |  y)$(w  |  y')|  d3y  d3y'.  (D.IO) 

The  double  integral  on  the  rhs  of  (D  10)  separates  into  a  product  of  two  integrals. 
Using  this  we  rewrite  (D.10)  in  the  form 

Sw2(iu;r,z)  <  ^p(y)(Sr.b.(uj;y))1/2|G(w,r  +  y,z)||$(u>|y)|d3y^  .  (D.ll) 

By  rearranging  the  integrand  in  (D.ll)  and  using  the  fact  that  M^lfi(y)  is  a  prob¬ 
ability  density  we  can  write 

Swj(w;r,z)  <  Ml  (^[(|5t.b.(a>; y)||G(w,r  +  y,z)|2|§(a’  |  y)|J)l/J])  ,  (D.12) 

where  E[  ]  is  the  expectation  with  respect  to  the  probability  density  M^"lp(y)- 
Applying  the  inequality  (D,4b)  to  (D.12)  results  in 

Sw2(u>;r,z)  <  Mb £[|5r.b.(u;; y )||G(u>,r  +  y, z)|2|$(w  |  y)|2].  (D.13) 

Using  (D.2b)  to  bound  the  rhs  of  (D.13)  we  obtain 

Swj(w;r,z)  <  Afb£[|5r  b.(u;; y)||G(w,r  +  y,z)|2]|$  (w  |  y0(u>))|2,  (D.14) 

Finally,  by  expressing  the  integral  on  the  rhs  of  (D.lb)  in  terms  of  the  expected 
value  operator  £[■),  w«  obtain 

Swt(w;r,z)  =  Afbf/i|5rb,(w;  y)i|G'(u/,r  +  y,z)|J].  (D.15) 

The  upper  bound  follows  immediately  by  using  (D.15)  in  (D.14). 
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Appendix  E 

The  model  equations  for  the  examples  of  Sect.  5 

The  breaking-wave  energy  spectrum  and  the  noise  intensity  spectrum  describe  the 
noise  observed  on  a  single  hydrophone  located  at  the  position  z.  In  the  simplified 
model  for  the  ideal  medium  (homogeneous,  semi-infinite),  the  breaking-wave  energy 
spectrum  is  given  by 


5w(o>;r,z)  =  MbSm.b.(w)ii(w,0):r(w1|z  -  rj),  (E.la) 


and  the  noise  intensity  spectrum  is  given  by 


S(w,  z)  -  AMbSm.b.(w)<?'(u>;.t,0,7).  (E.lb) 

The  transmission  loss  T(u>, |z  -  r|)  is  determined  by  (14c)  and  the  integrated  spec¬ 
trum  C-(im;  2,0,7)  is  determined  by  (15b)  with  the  hydrophone  separation  z'  set 
equal  to  zero.  The  mean  bubble  spectrum  Sm ,b.(a>)  and  the  radiation  pattern  R(u>,  9) 
depend  or.  the  mechanism  by  which  the  bubble  vibrations  are  excited.  In  this  ap¬ 
pendix  we  derive  the  equations  that  determine  the  mean  bubble  spectrum  and  the 
radiation  pattern  for  the  excitation  postulates  of  Sect.  4  and  the  assumption  that 
the  pressure  transfer  spectrum  H(u>\  y,  R),  the  bubble  generation  rate  p(y)  and  the 
bubble-radius  density  p{R  |  y)  are  independent  of  y  throughout  the  source  region. 

For  the  pressure-jump  excitation  the  radius-averaged  bubble  spectrum  Srh(u>;y) 
is  obtained  from  (8c)  with  the  single-bubble  spectrum  5b(i*,;y,J?)  obtained  from 
(Sd)  using  the  excitation  energy  spectrum  of  (18b).  Making  these  substitutions  and 
noting  that  Si?{u,y ,  R)  is  independent  of  R  results  in 

Sr.b.(w;y)=  E[\H(u>,R)\2}(pgd(y)/w)?,  (E.2) 

where  £?[•]  is  the  expectation  with  respect  to  p(fi).  The  mean  bubble  spectrum  is 
obtained  by  using  (E.2)  in  (8b).  The  result  is 

Sm.b  (w)  =  E[\H(w,R)\2}(pad{y)/^)iM^1p(y)d3y 

=  E\\a(*r,  Jt)|*](WM*  (  (d(y))JAfb~W)dV 

Jv 


or  equivalently 


5m.b.H  =  5r.b.(w;X>').  (E.3a) 

where  ST  b.(w;  D1)  is  the  radius- averaged  bubble  spectrum  evaluated  at  <f(y)  — -  D' , 
and 


(E.3b) 
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is  the  rms  source  depth. 

The  normalized  radius-averaged  spectrum  5rb .(<*»; y )  is  obtained  by  using  (E.2)  in 
(8f).  The  result  is 

■Sr.b.(w;y)  =  Wb-V(y)£[|WKi?)|2l(psd(y)M2/5mb.(u») 

or,  alternatively,  using  Eqs.  (E.3) 

•5r.b.(u»;y)  -  (E.4) 

In  the  simplified  model  the  source  region  has  a  horizontal  cross-section  A1  that  is 
independent  of  depth  so  that 

d(  y)  =  y",  (E.5a) 

and  for  the  special  case  where  /i(y)  is  constant 

My)  =  Mh/(A'D).  (E.5b) 

It  follows  from  (E.3b)  and  Eqs.  (E.5)  that  D‘  is  given  by 

'  D  \1/5 

(»")2dy"l  =  0.58D,  (E.6) 

and  from  (12c),  (E.4)  and  Eqs.  (E.5)  that  5r.b.(w;y")  is 

5,.b.(u>;y")=  D-V)*/(0')*.  (E.7) 

The  radiation  pattern  is  obtained  by  using  (E.7)  in  (14b).  The  resu’t  is 

R(u>,8)  =  D~l(D')~1  fD(y")2  (2sin((u;/c)co *(8)y"))*  dy" .  (E.8) 

Jo 

For  the  volume-rate  excitation  the  radius-averaged  bubble  spectrum  is  obtained  from 
(8c),  (8d)  and  (19b).  The  result  is 

5rb(w;y)=  E[|/f(u;;JR)|J/?3l(pl/')2.  (E.9) 

Since  the  rhs  of  (E.9)  is  independent  of  y,  the  mean  bubble  spectrum  is  also  given 
by  (E.9)  and  the  normalized  radius- averaged  spectrum  is  simply 

S».b.(w;y)  =  Mb“‘My)-  (E.10) 

For  the  simplified  model  with  /i(y)  constant  throughout  the  source  region,  it  follows 
from  (12c),  (E.10)  and  Eqs.  (E.5)  that  St.b.(w;  y")  is 

S.bV.y")  =  D-\  (Ell) 

and  from  (14b)  and  (E.ll)  that  the  radiation  pattern  is 

[D 

ft(w,0)=D-‘  /  (2 sin ((w/c) cos(0)y"))2  dy”.  (E.12) 

Jo 
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